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We propose steady-state electron transport based on coherent transfer by adiabatic passage 
(CTAP) in a linearly arranged triple quantum dot with leads attached to the outer dots. Its main 
feature is repeated steering of single electrons from the first dot to the last dot without relevant 
occupation of the middle dot. The coupling to leads enables a steady-state current, whose shot 
noise is significantly suppressed provided that the CTAP protocol performs properly. This repre- 
sents an indication for the direct transfer between spatially separated dots and, thus, may resolve 
the problem of finding experimental evidence for the non-occupation of the middle dot. 
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The possibility of finding particles in delocalized states 
is a most intriguing feature of quantum systems. It gives 
rise to phenomena such as tunneling which is intrinsically 
quantum mechanical and beyond everyday experience. A 
particular tunnel effect is the coherent transfer by adi- 
abatic passage (CTAP) [TJ of an electron in a linearly 
arranged triple quantum dot [3]. It is based on the 
control of the inter-dot tunnelings, such that an electron 
is guided from the first dot to the last dot without pop- 
ulating the middle dot. While being closely related to 
stimulated Raman adiabatic passage [U [5] , CTAP pos- 
sesses the appealing feature of dealing with spatially sep- 
arated states. 

The non-occupation of the middle dot can be the- 
oretically predicted straightforwardly by solving the 
Schrodinger equation or, in the presence of decoherence, 
an appropriate quantum master equation. By contrast, 
its direct experimental proof is less obvious, because such 
position measurement unavoidably acts back on the po- 
sition. Thus, such measurement would influence the ef- 
fect that it should substantiate. For example, coupling 
the middle dot to a charge detector [5J [7] impacts signifi- 
cantly on the dot populations [8] , since it causes decoher- 
ence to which CTAP is quite sensitive ,.9 . Consequently, 
as a matter of principle, direct observation of the dot 
populations is of limited use and, thus, one may have to 
resort to an indirect verification. 

In this work, we propose a generalization of CTAP that 
circumvents this difficulty and predicts an experimental 
fingerprint for the non-occupation of the middle dot. It 
is based on a modified setup with electron reservoirs at- 
tached to the first and the last quantum dot (see Fig. [IJ , 
such that a current may flow. We demonstrate that 
the suggested steady-state CTAP protocol establishes a 
rather regular electron transport with suppressed shot 
noise indicating its proper course. In turn, if the mid- 
dle dot is significantly occupied, transport becomes more 
noisy. Our analysis of the transport process is based on 
a numerical propagation of a time-dependent quantum 
master equation, which needs to be augmented in order 
to compute current noise jTUl [TTJ . 



Model and master equation. — The setup sketched in 
Fig.[l]is described by the Hamiltonian H(t) = -£/tqd(£) + 
-Hdot-ieads + H r, + Hr, where (in units with h = 1) 

H TQD (t) = n 12 {t)c\c 2 + n 23 (t)clc 3 + H.c. (1) 

refers to the triple quantum dot with vanishing onsite 
energies. We assume strong Coulomb repulsion such that 
the triple dot can be occupied with at most one electron. 
Hence the electron spin has no relevant influence and will 
be ignored. By applying appropriate gate voltages to the 
inter-dot barriers, the tunnel matrix elements between 
neighboring dots, Qi 2 (t) and ^23^), are controlled such 
that they assume the sequence of double Gauss pulses 



12/23 



(*) 



OO 

E 

n=0 



(t± At/2 -nT -t ) 2 
2cr 2 



(2) 

with amplitude f2 ma x, width 2a, delay time At, and rep- 
etition time T as is depicted in Fig. [2]ja). This implies 
the time periodicity H (t) = H(t + T) for t > to with 
some time-offest to- 

Ideally, one would like to modify the control parame- 
ters Clij adiabatically slowly to ensure ideal accomplish- 
ment of the protocol. This however is not possible, be- 
cause a lower limit to the CTAP duration is set by de- 
coherence, which may stem from substrate phonons, a 
detector, or as in the present case the coupling to leads. 
For the transient dynamics during a total propagation 
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FIG. 1. Triple quantum dot in linear arrangement, where the 
first and the last dot are tunnel coupled to electron reservoirs 
with a large voltage bias applied. The CTAP protocol can be 
implemented by proper time-dependent gate voltages at the 
inter-dot barriers, such that the corresponding tunnel matrix 
elements f2i2(£) and f223(t) obey Eq. 
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FIG. 2. (Color online) (a) Time profile of the inter-dot tunnel 
matrix elements forming repeated CTAP pulses. (b,c) Time 
evolution of the dot occupation probabilities for T = 0.5/f2 max 
and driving periods (b) T = 2/f2 max and (c) T = 2.9/f2 ma x- 



time i max considered in Ref . [I] , decoherence is tolerable 
for the pulse parameters At = 2cr = T/4. Here we use 
T = frmax as pulse repetition time and focus on the steady 
state. An important feature of the CTAP double pulse 
is its "counter-intuitive" order pQ, which means that the 
coupling 2 3 is active before r2 12 . This requirement can 
be understood upon noticing that (^23, 0, — f2i2) T is an 
eigenstate of the single-particle Hamiltonian in the basis 
{|1), 1 2), 1 3)}. Between the maxima of the double peak, 
this eigenstate changes from (1,0, 0) T to (0, 0, — 1) T . 
Thus, the adiabatic time evolution of an electron start- 
ing in dot 1 is as deserved, namely it undergoes a direct 
transition to dot 3 without populating dot 2. 

The leads are modeled as free electron gases with the 
Hamiltonian i/i ca d s = Ylt q e q c iq c tq> where c\ q creates an 
electron in state q of lead £ = L,R. The coupling of 
dot 1 to the left lead is given by the tunnel Hamilto- 
nian H l — J2 q VLqC^qCi + H.c. and Hr accordingly. The 
lead coupling is fully determined by its spectral density 
rt(e) = 27r^ 9 \VL q \ 2 8(e — e q ), which we assume within a 
wide-band limit energy independent and equal for both 
leads, T L (e) = T R (e) = T. 

If all relevant dot energies lie within the voltage win- 
dow of the two leads, electron transport becomes unidi- 
rectional. Then one can derive for the reduced triple dot 
density operator p within a Bloch-Redficld approach the 
master equation |12) 

p= -i[H TQD (t),p]+TV(c\)p + rV(c 3 )p = C(t)p (3) 



with the Lindblad form T>(x)p = xpx 
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The dissipative superoperators 23(c{) and 'D[c 3 ) describe 
incoherent tunneling from the left lead (source) to dot 1 
and from dot 3 to the right lead (drain), respectively. 

Computing current noise by numerical propagation. — 
The fluctuations of the current can be characterized by 
the current auto correlation function S{uj) which is typi- 
cally measured at low frequencies. Therefore, we focus on 
the zero-frequency limit in which the correlation function 
relates to the variance of the transported charge accord- 
ing to S = 5(0) = e 2 lim t ^ 00 (d/dt)(AN 2 (t)) [Ml E], 
where e is the elementary charge. The variance can be 
computed with the help of a density operator R{x, t) that 
is augmented by counting variable x an d possesses the 
limit iJ(x->0,t) = p(t). Its trace trR(x,t) = (e lxN ) 
is the moment generating function for the number of 
transported electrons N, i.e., (N k ) — (d/dix) k tr R\ x —q. 
For unidirectional transport, R obeys the master equa- 
tion R(x,t) = CR(x,t) + (e l * - l)JR( X ,t) [TS], where 
JR = Tc^Rcl is the jump operator contained in the 
Lindblad form. Exploiting the 2ir periodicity in Xi 
one frequently continues with a Fourier transformation 
which yields the so-called JV-resolved master equation 
|12) . Then however, a major technical problem arises, be- 
cause one has to consider many additional degrees of free- 
dom, such that a numerical solution becomes rather ex- 
pensive. A more economic approach has been developed 
in Refs. [TU1HI]- By Taylor expansion up to second order, 
-R(x) — P + ixRi ~ |x 2 -^2) one obtains the equations of 
motion Ri = CRi + Jp and R 2 = CR 2 + 2JR 1 + J p, 
which we solve by numerical propagation. Then we use 
the fact that tiR(x) is the moment generating function, 
so that (d/dt)(N k ) = tr Rk, by which we finally obtain 
the current I = e lim t _i. 00 tr R\ and the zero-frequency 
noise S = el + lim^oo 2e tr(ej'i?i — IRi). For transport 
in periodically time-dependent conductors, these quanti- 
ties must be considered in the average over the driving 
period [IB] . 

The relative noise strength is characterized by the Fano 
factor F = S7e|/| which provides information about the 
transport mechanism [T3]. For a Poisson process it has 
value one, while it becomes smaller for more regular 
transport. In the present context, we will speak of "sig- 
nificant shot noise suppression" if the Fano factor lies 
clearly below 1/2. 

Time evolution. — In order to acquire insight into the 
system dynamics, we time integrate master equation ^ 
by a Runge-Kutta method. We start with an empty triple 
quantum dot just before the first double pulse sets in 
and obtain the time evolution of the dot occupations de- 
picted in Fig. [2pD). At the initial stage, t < 0.1T, dot 1 
becomes populated by an electron that tunnels in from 
the left lead. Between the peaks of the double pulse 
[cf. panel (a)], the electron essentially undergoes a di- 
rect transfer to dot 3. As non-adiabatic correction, dot 2 
acquires some occupation, which however stays clearly 



3 



(a) 3 



a 2 
a 



(b) 3 




_ 0.5 




[1/Qn 



FIG. 3. (Color online) (a) Time averaged population of dot 2, 
P22, and (b) Fano factor F = S/e\I\ in the steady-state limit 
as function the pulse repetition time T and the dot-lead cou- 
pling r. For graphical reasons, the data are clipped at the 
values 0.5 and 1, respectively. 



below the occupation of the other dots. After the dou- 
ble pulse, the electron tunnels further to the right lead, 
and a next electron may enter at dot 1. The fast oscil- 
lation of P22 when inter-dot tunneling is active can be 
identified as tunnel oscillation between the middle dot 
and one of its neighbors, because its frequency is roughly 
^max/27r [T7] ■ For certain values of the driving period 
T [panel (c)], these oscillations become larger and may 
interfere constructively, so that dot 2 acquires a larger 
occupation. 

If all tunnel matrix elements were constant in time, 
the population of all three dots were about equal (not 
shown). Thus, our transport process must be dominated 
by CTAP if the occupation of the middle dot is consider- 
ably smaller than the occupation of the other two dots. 
In the ideal case, the steady state will be such that with 
each double pulse, one electron is transported. This gives 
rise to a regular current, so that shot noise is suppressed. 
In the following, we analyze this current and will find 
that the proper course of the CTAP protocol correlates 
with small noise manifest in the Fano factor. 

Steady state. — For computing the current noise, we 
need to solve simultaneously master equation ^ and the 
equation of motion for Rx, as discussed above. Since we 
are interested in the long-time behavior, we eliminate 
transient effects by considering the results only after all 
elements of the density matrix have changed by less than 
0.25% with respect to the previous period. Convergence 
is typically reached after five to ten double pulses. 

Figure [3] provides an overview of the global behav- 
ior as function of the driving period T and the dot-lead 
coupling T. For the time- averaged occupation of dot 2 
shown in panel (a), we distinguish three regions. First, 
for T < l-5/O max in which the driving is too fast, such 
that non-adiabatic effects play a strong role, and we can- 



not expect CTAP to work properly. Accordingly, the 
average population P22 may take any value. The corre- 
sponding Fano factor [panel (b)] is about 1/2 or larger, 
i.e., there is no relevant shot noise suppression. Notice 
the similarity of the dot occupation with the one ob- 
served in Ref. PQ, despite that there transient dynamics 
with phenomenological decoherence was considered. 

For larger driving periods, we find a broad region in 
which the population of the middle dot is rather low, 
unless the dot-lead coupling T becomes too large. In 
particular, we find islands in which both the average 
and the maximum of the mean occupation is as small 
as P22 ~ P22 ~ 0.1. The plot for the Fano factor exhibits 
a similar structure, where both F and P22 assume small 
values for roughly the same parameters. This provides a 
first hint that our generalized CTAP gives rise to regu- 
lar, low-noise electron transport, in particular when T is 
of the order f2 max - Interestingly, for certain values of T, 
dot 2 becomes occupied and the Fano factor increases ac- 
cordingly. This can be attributed to the mentioned tun- 
nel oscillations, because the distance between two peaks 
increases linearly with the amplitude f2 max . 

For a more quantitative analysis, we depict in Fig. [4] 
the time-averaged steady-state occupation p 2 2 and the 
Fano factor for two relatively small values of the dot- 
lead coupling r. This demonstrates that the population 
of dot 2 can be as low as 0.1, while it becomes roughly 1/3 
in the intermediate region. Accordingly, the Fano factor 
changes from a significantly sub-Poissonian value toFw 
1. The dc current [inset of Fig.^a)], by contrast, exhibits 
only minor signatures of CTAP. This underlines that it 
is indeed the noise which provides information about the 
nature of the transport process. Moreover, it reveals that 
for the relatively small T considered, on average only each 
second double pulse transports an electron. This value 
can be improved by increasing T which, however, requires 
a smaller dot-lead tunneling T. 

We quantify the relation between the population P22 
and the Fano factor F by the correlation coefficient 
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between these quantities for T in the relevant interval 
[2/fi max ,4/fi max ]. The numerator denotes the covari- 
ance, while the denominator contains the standard de- 
viations, such that r = 1 for ideal correlation. The inset 
of Fig. |4|b) shows that for T < !] max , the correlation 
is rather high, typically 0.75 or larger. Moreover, the 
difference between the maximum and the minimum of 
the Fano factor, F max — F min , is roughly 0.5, i.e., suffi- 
ciently large to be discriminated in an experiment [18] . 
For larger values of T, the Fano factor and the popula- 
tion become significantly anti-correlated. There however, 
-Fmax — -Fmin is quite small and may be below the exper- 
imental resolution. Therefore we do not further discuss 
the regime T > 1.5£l max . 
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FIG. 4. Fano factor F — S/e\I\ and time-averaged popu- 
lation P22 for the dot-lead couplings (a) V = f2 m ax and (b) 
r = 0.5f2 m ax- The horizontal dotted line marks the value 0.1 
which is the minimal reminiscent population of dot 2. In- 
sets: (a) dc current corresponding to the main plot in units of 
electron per pulse, (b) Correlation coefficient p} of the Fano 
factor and P22 in the range [2/f2 max , 4/f2 max ] (solid line) and 
the corresponding minimum and maximum of the Fano factor 
(dashed). 



In a possible experimental realization of steady-state 
CTAP, the oscillatory dependence of the Fano factor as 
function of the pulse distance T may serve as indica- 
tor. Most significant results are expected for T < fl max . 
Thus, one may use a triple quantum dot whose four bar- 
riers have similar transparency. Negative gate voltages 
applied to the inner barriers allow one to control the 
tunnel matrix elements. Since the average current con- 
sists of roughly half electron per double pulse, a current 
/ m 10 pA requires a repetition rate of 100 MHz [19]. The 
maximal tunnel rates should be considerably larger than 
that, say, by a factor 10, which means that T and O max 
should be of order 10 fieV. 

Conclusions. — We have proposed a CTAP protocol for 
a triple quantum dot in transport configuration, i.e., 
when the first dot and the last dot are coupled to electron 
source and drain. It induces a steady-state transport in 
which electrons tunnel non-locally from the first to the 
last dot, practically without populating the middle dot. 
The noise properties of the resulting current depend on 
the size of the reminiscent population. The difference is 



manifested in the Fano factor which, thus, represents a 
fingerprint of the proper course of the protocol. Such in- 
direct evidence for CTAP is particularly useful, because 
any direct measurement of the dot occupations would 
act back on the system and, thus, unavoidably influence 
the effect to be substantiated. A particular feature of 
steady-state CTAP is that it fails for certain distances 
between pulse pairs. Thus, upon changing this distance, 
one obtains a sequence of sizable shot noise suppression 
and shot noise enhancement, which should be measurable 
by standard techniques. Thus, our proposal may initiate 
the experimental realization of this intriguing non-local 
transport phenomenon. 
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